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Abstract 

We investigate the effects of a squeezed pump on the quantum properties 
and conversion efficiency of the hght produced in single-pass second har- 
monic generation. Using stochastic integration of the two-mode equations 
of motion in the positive-P representation, we find that larger violations of 
continuous- variable harmonic entanglement criteria are available for lesser ef- 
fective interaction strengths than with a coherent pump. This enhancement 
of the quantum properties also applies to violations of the Reid-Drummond 
inequalities used to demonstrate a harmonic version of the Einstein-Podolsky- 
Rosen paradox. We find that the conversion efficiency is largely unchanged 
except for very low pump intensities and high levels of squeezing. 
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1. Introduction 



One of the simplest non-linear optical processes is travelling wave second 
harmonic generation lj|, in which a nonlinear x^"^^ crystal is pumped with 
an electromagnetic field at one frequency and produces a second harmonic 
field at twice this frequency. A comprehensive classical treatment of this pro- 
cess was first given by Armstrong et al. The quantum properties of the 
output fields were first calculated using a method of linearisation about the 

n 

classical solutions [3[, even though it had long been known that these were 
not accurate for arbitrary interaction strengths 4] • The approach of Walls 
and Tindle, using matrix equations for the number state coefficients, was 
necessarily limited to small photon numbers. The later development of the 
positive-P representation sj], which maps the quantum evolution equations 
onto stochastic differential equations in a doubled phase space, allowed for 
the treatment of much larger photon numbers. For this system, this was first 
taken advantage of by Olsen et al. 6| to treat the two-mode model, find- 
ing that full conversion to the second harmonic did not occur, but that the 
fundamental field experienced a revival inside the crystal. This approach en- 
abled a calculation of the quantum properties, such as quadrature squeezing, 
of the output fields without relying on any assumptions about the mean-field 
solutions. The positive-P representation was also used to calculate the QND 
(Quantum Non-Demolition) properties of the system j?] and also compared 
with results found using the semi-classical method theory of stochastic elec- 



trodynamics [8[. The quantum correlations between the two fields were also 
calculated {q], using correlations which later became famous as the Duan- 



Simon criteria for two-mode continuous- variable entanglement 10|, lUl]- The 
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entanglement between the fundamental and harmonic fields was later named 



harmonic entanglement by Grosse et al. 



12| . who calculated it for a system 



which could operate in both the up and down-conversion regimes. The Reid 



correlations [IS] for Einstein-Podolsky-Rosen entanglement [IJ] between the 
two modes have also been calculated previously 'l5|. In the spirit of Grosse 
et ai, we shall name this harmonic steering. 

The idea that the conversion efficiency in nonlinear optical processes couk 
be a function of the quantum statistics of the inputs was raised by Shen jl6| 
in 1967. Shen showed that the conversion efficiency in the two-mode model 
of second harmo nic g eneration would depend on the second-order correlation 
function, g^^\0) (l7|, of the pump, predicting that light with chaotic statis- 
tics would initially convert twice as efficiently as a coherent pump. This was 
later verified by stochastic integration, where the differences in efficiency be- 
tween these pumps were explicitly calculated [l8| without making the small 
interaction approximations required in Shen's approach. The development of 
algorithms to model different quantum states in the positive-P and Wigner 



representations 



191 ] allowed for the investigation of the effects of these in 



the process of atom-molecule conversion in trapped Bose-Elnstein conden- 
sates 
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2l| . showing that the initial quantum statistics would also affect 



this process. Because squeezed states with low average photon number have 
(7*^^^(0) > 1, we decided to investigate the effects of input squeezing in this 
system. In this work, we use the method given in Ref. 19| to model the in- 
put to our crystal as squeezed states with different degrees of amplitude and 
phase quadrature squeezing, and investigate the outputs in terms of conver- 
sion efficiency, single- mode squeezing, two- mode entanglement, and the EPR 
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paradox as expressed by the Reid criteria. 

2. Squeezed states 

Before we enter into the calculations for the problem, we will review the 
definition of squeezed states and the quadratures as we define them in this 
work. This is necessary as there are several different conventions found in the 
literature, which can lead to confusion if the quantities used are not explicitly 
defined. We begin with the bosonic annihilation operators, d for photons in 
the fundamental mode at frequency uj, from which we define the quadrature 
operators 

Xa^a + a), Ya = -i{a-a'), (1) 

with similar definitions for the harmonic mode at 2uj, using the operator b. 
With this definition of the quadratures, the vacuum or coherent state level 
of the variances is 

V{X,) = V{Y,) = 1, (2) 

where j — a,b. 

A squeezed state is then a state of the electromagnetic field that has a 
variance of less than one in one of the quadratures, at the expense of increased 
fluctuations in the orthogonal quadrature. We note here that the quadratures 
can be defined at any angles, but that zero and 7r/2 are sufficient for our 
purposes here. A coherently displaced squeezed state is written as |Q;,re'^), 
where the c-number a represents the coherent displacement and r is known as 
the squeezing parameter, with the angle of the squeezed quadrature. This 
state has a mean intensity, Na = + sinh^r and quadrature variances for 
= 0, V{Xa) = e-^"^ and V{Ya) = e^^ For = n/2, the Y quadrature is 
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squeezed and X is anti-squeezed. The formal squeezing operator is S{e) 



exp(l/2e*a2 — l/2ea^2) 22|, where e = re^***^. As this operator creates and 
annihilates photons in pairs, we might expect a squeezed field to lead to 
increased second harmonic conversion efficiency, which we will investigate in 
what follows. 



The second-order correlation function is 



(3) 



which can be expressed in terms of the number variance as 



We can immediately see that if we have a super-Poissonian field with V{N) > 
(a^a), g^'^\0) will be greater than unity and we will find photon bunching. 
On the other hand, V{N) < {p)a) will give an anti-bunched field. 
The explicit expressions for the number variances are found as [l| 

V"(A^) = |ap exp(— 2r) + 2 sinh^ r cosh^ r, X squeezed 

V{N) = lap exp(2r) + 2 sinh^ r cosh^ r, Y squeezed, (5) 

from which we see that a field squeezed in the X quadrature will tend to 
exhibit anti-bunching while one squeezed in the Y quadrature will tend to 
exhibit bunching. However, calculations show that even for a coherent ampli- 
tude of a = 10^, the values of g^'^\Q) remain close to unity. As the coherent 
amplitude increases, they diverge even less from unity, so that we would not 
expect a squeezed pump to result in significantly enhanced conversion effi- 
ciency except for very weak pumps. We will investigate the effects on the 
quantum properties of the output ffelds below, using stochastic integration. 



3. Hamiltonian and Equations of Motion 

In this work we will use a simplified description of travelling wave second 
harmonic generation which does not treat dispersion within the x*^^^ medium 
and treats the fields as plane waves at fixed frequencies, u and 2u. Whilst 
this approach would not be highly accurate if we sought to model an actual 
experiment, it is sufficient to show the effects we are looking for as functions 
of the degree of squeezing of the pump beam. The interaction Hamiltonian 
is written as 



H = ih- 

2 



(6) 



where k is the effective nonlinearity. Following the standard procedures |23| . 
we map this Hamiltonian onto a set of four stochastic differential equations 
for the evolution of the c-number variables of the positive-P phase-space 
representation as the fields traverse the nonlinear medium. Proceeding via 
the von Neumann and Fokker-Planck equations, we find 

da 

dz 
da~^ 



dz 
dp 



Ka/5+ + a/ Kl3+rj2{z), 



— a 
dz 2 

dP+ K +2 
= ry^ 

dz 2 



(7) 



where the rij{z) are real Gaussian noise terms with the properties rij{z) = 



and ■r]j{z)rik{z') = 6jkS{z — z'). The independence of the two noise terms 
means that a and are not complex conjugate except on average, which also 
holds for P and P'^. It is this property which allows us to integrate what are 
equivalent on average to equations of motion for non-commuting operators. 
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Averages over a large number of trajectories of the system of equations above 
allow us to find normally-ordered operator expectation values, with 

ama+n ^ (at"a"^), (8) 

and similarly for 6, and /3, (3^ . In practice we integrate these equations in 
Matlab and average over at least 10^ trajectories. 



4. Squeezing for Different Input States 

We have investigated the outputs for pump fields which are squeezed in 
either the X (amplitude) or Y (phase) quadratures, for squeezing parame- 
ters, r = 0, .5, 1. The results are shown in the graphs as a function of the 



parametrised interaction length, C, = Kzy (|a(0)p + sinh^r)/2. In all our 
calculations in this section, we use a(0) = 10'^ and k = 10~^. For these input 
parameters, we find no significant difference in the conversion efficiencies, 
consistent with the fact that ^'^^■'(O) stays very close to unity. We did find 
that for even smaller a(0), the efficiency did improve more with increasing Y 
squeezing, but as we are interested in the quantum properties of reasonably 
intense outputs, we have not shown these results here. 

The first quantum correlations we investigate are the quadrature vari- 
ances, since bright squeezing can have applications in quantum communi- 



cation and and teleportation 



24 



25| . Beginning with a coherent pump. 



the equations of motion derived from our Hamiltonian predict almost per- 
fect squeezing in the fundamental after a certain interaction length, with a 
maximum of 50% squeezing in the harmonic, with the X quadrature being 
squeezed in both cases |6i]. After obtaining a maximum of squeezing, both 
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Figure 1: V{Xa) for different input states. On the left Irand side of the plot, starting 
from the bottom, we show the values for X-squeezed with r = 1 (dashed line), r — .5 
(dash-dotted line), r = (solid line), then y-squeezed with r — .5 (dotted line) and r = 1 
(dashed line with crosses). We note here that the method used to perform the integration 
does not allow for an easy determination of the sampling errors, but when we run the 
program several times and compare the results, we do not find significant differences. For 
the quantities shown here, the inferred errors would be of the size of the line widths. The 
quantities plotted in this and all subsequent plots are dimensionless. 

become antisqueezed as the fundamental mode revives inside the nonhnear 
medium. 

In Fig. [T]we show the pump quadrature variances, V{Xa), for different 
input statistics, with the sohd hne beginning at one being the prediction for a 
coherent pump. We see that a pump squeezed in the X quadrature actually 
results in output fields which are squeezed to begin with, as expected, but 
that they lose some of this squeezing until they reach the same level as that 
for a coherent pump, which happens just beyond C = 1- After this they 
become increasingly squeezed, staying below the level for the coherent pump 
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at least until the squeezing begins to diminish, at C ~ 4.5. For a pump 
with Y squeezing, we see that the output fields begin as antisqueezed, with 
squeezing developing until they are equal to the other fields shown around 
( ^ 1.2. The variance of the field resulting from the pump with r = 0.5 
continues to decrease, but stays above the level for the coherent input until 
C ^ 5. Beyond the length where the variances are equal, it begins to increase 
for r = 1, but then decreases again until the variances of the five different 
fields are indistinguishable at around ( ^ 5. This behaviour shows that 
the system has some memory of the input even when the fundamental field 
has reached its minimum and begun to revive, which happens in the region 
4 < C ^ 6 for the parameters we use here j6|. The quadrature always 
became antisqueezed, even when it began with V{Ya) < 1 as a result of a 
Y squeezed pump. The output quadrature was never more squeezed than 
the input quadrature, with the variances increasing monotonically over the 
range that we investigated. 

When we look at the variances in the harmonic, we find interesting be- 
haviour, which is qualitatively different depending on which pump quadra- 
ture is squeezed. The effects on V{Xb) of squeezed pumps are shown in Fig.[2l 
For a coherent pump we see that the squeezing in this output quadrature 
reaches a maximum at V{Xb) ~ 0.5, with Xa-squeezed pumps resulting in 
a significantly greater degree of squeezing, while Fa-squeezing in the pump 
results in antisqueezing in the Xf, output. However, when we look at the 
behaviour of the Yb quadrature, we see that a Ya (phase) squeezed pump 
leads to squeezing in this Yb quadrature, as well as in Xa- This is in stark 
contrast to an amplitude quadrature squeezed pump, which gives squeezing 
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Figure 2: V{Xb) for different input states. On the left hand side of the plot, starting 
from the bottom, we show the values for X-squeezed with r = 1 (dashed line), r — .5 
(dash-dotted line), r = (solid line), then F-squeezed with r — .5 (dotted line) and r — 1 
(uppermost dashed line with crosses). 

in both the output amphtude quadratures. This effect may be of interest ex- 
perimentally for controlling the harmonic quadrature in which squeezing is 
available, as this changes from amplitude to phase quadrature as the squeezed 
pump quadrature changes, while leaving the fundamental squeezing in the 
Xa quadrature. 



5. Entanglement and EPR for Different Input States 



Squeezing in and of itself is a single-mode phenomenon, and even though 
it has several potential applications, the field becomes more interesting when 
we consider multi-mode correlations, defined by combining the natural modes 
of the problem in various manners. For the present system, this allows us to 



consider continuous-variable bipartite entanglement 
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ll| and demonstra- 
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Figure 3: V{Yi,) for pump fields initially squeezed in the Y quadrature. V{Ya) begins at 
the initial pump values and rapidly becomes anti-squeezed. With an X-squeezed pump, 
V{Yb) always shows noise above the coherent state level. 



tions of the Einstein-Podolsky-Rosen paradox 



more recently labelled 



steering [26.]. Asymmetric Gaussian steering has previously been theoreti 
cally analysed using an unbalanced nonlinear coupler ^\ and experimen- 



tally produced using parametric downconversion 



28j. The phenomenon of 



entanglement between the fundamental and harmonic fields has been pre 



dieted 9| and subsequently labelled harmonic entanglement 12|. Following 
in this spirit, we will use the name harmonic EPR for the EPR correlations 
we calculate between the fundamental and harmonic fields. 

We will use the Duan-Simon inequalities in their simplest forms, where 



ViXa±X,) + V{YaTYt,)<^ 



(9) 



guarantees that field a and field b are entangled. We note here that a value of 
this inequality greater than or equal to 4 does not prove that entanglement 
is not present, as we have not taken advantage of the optimisation avail- 



Figure 4: The Duan-Simon correlations with the pump squeezed in the X quadrature. 
SoHd hne is r = 0, dash-dotted Hne is r = 0.5 and dashed line is r = 1. The straight line 
gives the value expected for uncorrelated coherent states. 

able in the original definition. This is not important for our purposes here, 
because we will later examine the EPR/steering correlations, which are of 
more practical significance for quantum information tasks. It is easily seen 
by inspection that two uncorrelated coherent states will give a value of 4 for 
the expression in Eq. [91 

In Fig. m we show the results obtained when the pump is X-quadrature 
squeezed, with the solid line being for an initial coherent state with r = 0, 
the dash-dotted line representing r = 0.5, and the dashed line representing 
r = 1.0. We readily see that increasing the input squeezing increases the 
degree of violation of the inequality, while also causing the violation to begin 
for smaller interactions. In this case, the violation is found by measuring the 
combined quadrature variance V{Xa — Xb) + V(Ya + Yb). In Fig. |5]we show 
the effects of input squeezing in the F-quadrature. We immediately see that 
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Figure 5: The Duan-Simon correlations with the pump squeezed in the ^-quadrature. 
Sohd hne is r = 0, dash-dotted Hne is r = 0.5 and dashed line is r = 1. The straight line 
gives the value expected for uncorrelated coherent states. 



in this case, the violation is found by measuring the combined quadrature 
variance V{Xa + ^b) + V{Ya — Yfy), which is not violated for a coherent 
pump. The lines for different squeezing values are as in Fig. |H Although not 
investigated here, it seems reasonable to infer that the violations will be found 
at different quadrature angles of the expression in Eq. |9] as a continuously 
changing and linear function of the angle of input squeezing. This tunability 
could possibly be of use in quantum information applications. 

In order to investigate harmonic EPR/steering we have extended the 
approach of Reid ^| to use two-frequency inferred inequalities as previously 
defined by Olsen 15|. The basis of the EPR paradox is that, by inferring the 
values of one quadrature by using measurements on a conjugate quadrature, 
and then reversing the operation, we can obtain inferred measurements with a 
product of uncertainties which seemingly violates the Heisenberg Uncertainty 
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Principle. Of course, as has been exhaustively discussed in the literature, 
there is no actual violation due to the fact that the quadrature values have no 
defined reality before the act of measurement. To see the paradox in action, 
we define inferred values of quadrature X^, Ya which may be obtained via 
measurements of quadrature Xb, Yj,, and vice versa. 



v{t) - 

V{Xt) - 



V{Xa,Xb 



1 2 



ViXb) 



ViYb) 

V{Xa,Xb 



ViXa) 
V{Ya,Yb) 



ViYa] 



(10) 



where V{A,B) = (AB) — {A){B). We see a demonstration of the paradox 
whenever one of the following holds 



V'^f{Xa)V{Ya) < 1, 

y*"^(Xfe)F*"^(r5) < 1. 



(11) 



We note here that this is different to the standard definition, where quadra- 
ture measurements at the same frequencies are used. In our case, we are using 
a quadrature measurement at the fundamental (harmonic) frequency to in- 
fer values at the harmonic (fundamental) frequency. Since the Heisenberg 
Uncertainty Principle holds for quadrature variances at the same frequen- 
cies, satisfaction of the inequalities of Eq. [11] provides a demonstration of the 
paradox. A satisfaction of one of the inequalities but not the other would 
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Figure 6: The harmonic EPR correlation V^"^^ {Xa)V'^^^ (Va) for an input pump squeezed 
in the X-quadrature. We note that the samphng error has increased in the region where 
the denominators obtain their minima, for the same numbers of trajectories as in previous 
figures. 

provide an example of asymmetric harmonic steering, but we did not find 
this here, with both the correlations of Eq. [11] returning identical values. 

In Fig. |6]we show the results for V^"'^ {Xa^V^"^^ (Ya) when we squeeze the 
input field in the X-quadrature. The continuous line is for a coherent pump, 
while the dash-dotted line is for a pump with r = 0.5 and the dashed line 
is for r = 1. Fig. [7] is the same except for the fact that the pump is now 
squeezed in the y-quadrature. We see that squeezing the pump in the X- 
quadrature gives a smaller maximum degree of the EPR correlation, although 
the values for smaller interaction strengths are increased. This may be of 
some benefit experimentally both because this is the region where we would 
expect our simulations to be more realistic and because they may be obtained 
with a shorter nonlinear crystal or less intense pump. In Fig. [7] we see that 
Y squeezing gives stronger steering at both shorter and longer interaction 
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Figure 7: The harmonic EPR correlation V^^^ {Xb)V"^^ (Yb) for an input pump squeezed 
in the K-quadrature. 

strengths, with only a small region where a coherent pump is superior. The 
maximum of steering also increases as the pump is squeezed more. This 
suggests that, if we want increased harmonic steering from this system, a Y- 
squeezed pump is superior apart from a small region. Since for a given crystal, 
the pump intensity can be tuned to stay away from this region, our results 
show clearly that pumping with a field squeezed in the F-quadrature is the 
best option. We note that the sampling errors increase in the regions where 
the denominators of the expressions are around their respective minima, but 
we have kept the same number of trajectories because the scaling of the errors 
changings with the square root of this number, meaning that we would need 
10^ to see one tenth of the error. The time required would be prohibitive. 
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6. Conclusions and Discussion 

We have used a simple model of single-pass travelling wave second har- 
monic generation with a squeezed pump to invetsigate the effects of the 
squeezing on the quantum correlations of the outputs. We found that for the 
reasonably intense pumping that we have used in our model, the conversion 
efficiency is basically unchanged by pump squeezing, with the second order 
correlation function of the input only rising marginally above that of a coher- 
ent state. When we investigate the quadrature variances of the output fields, 
we find that a squeezed pump can increase the squeezing available in the 
harmonic field drastically, while having less of an effect on the fundamental. 

The results are interesting when we look at the correlations between the 
fields at the two frequencies, in terms of both two-mode entanglement and 
EPR/steering. We find that use of a squeezed pump can increase the quan- 
tum correlations in the outputs, as well as making them available at lower 
interaction strengths. This couki lead to greater flexibility in the use of 
second harmonic generation, possibly the simplest nonlinear optical process, 
as a valuable tool in quantum information and communication science and 
technology. 
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